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(Dated: February 7, 2008) 

We measure the branching fraction for the flavor-changing neutral-current process B — > X s £ + £~ 
with a sample of 89 x 10 6 T(45) -> BB events recorded with the BABAR detector at the PEP-II 
e + e~ storage ring. The final state is reconstructed from e + e~ or fi~ pairs and a hadronic system 
X a consisting of one K or and up to two pions, with at most one We observe a signal 
of 40 ± 10(stat) ± 2(syst) events and extract the inclusive branching fraction B(B — > X s £ + £~) = 
(5.6 ± 1.5(stat) ± 0.6(exp syst) ± l.l(model syst)) x 10~ 6 for m t+l - > 0.2 GeV/c 2 . 

PACS numbers: 13.20.He, 12.15.Ji, 11.30.Er 



The rare decay B — > X s £ + £~ , which proceeds through 
the b — > s£ + £~ transition, is interesting because the study 
of its rate and charge asymmetry could lead to indi- 
rect observation of physics beyond the Standard Model 
(SM). This transition is forbidden at lowest order in the 
SM but is allowed at higher order via electroweak pen- 
guin and Vl^-box diagrams. This implies that non-SM 
physics in these loops would contribute at the same or- 
der as the SM [ij, |2| . Recent SM calculations of the in- 
clusive branching fractions predict B(B — ► X s e + e _ ) = 
(6.9 ± 1.0) x 1CT 6 [(4.2 ± 0.7) x 10~ 6 for m e+e - > 0.2 
GeV/c 2 ] and B(B -> X^+fj,-) = (4.2±0.7) x 10~ 6 00. 
Although the branching fraction measurement for inclu- 
sive decays is more challenging than for exclusive decays, 
it is motivated by smaller theoretical uncertainties. Ex- 
clusive B — > K^£ + £~ (£ = ejn) decays have been ob- 
served by Belle and BABAR 0, • Belle has also reported 
a measurement of the inclusive B — > X s £ + £~ branching 
fraction [||. 

The data sample used in this analysis was collected 
with the BABAR detector Q at the PEP-II asymmetric- 
energy e + e~ storage ring at the Stanford Linear Acceler- 
ator Center. The sample consists of 89 x 10 6 BB events 
recorded at the T(AS), corresponding to an integrated 
luminosity of 81.9 fb . 

We study the B — > X s £ + £~ process by fully recon- 
structing a subset of all possible final states. The re- 
constructed hadronic systems X s consist of one K ± or 
and up to two pions, with at most one tt°. This 
approach allows approximately half of the total inclu- 
sive rate to be reconstructed. If the fraction of modes 
containing a K L is assumed equal to that containing 
a K®, the missing states represent ~30% of the total 
rate. To compute the inclusive branching fraction, we 
account for missing modes and selection efficiencies us- 
ing a B X s £ + £~ decay model constructed as follows. 
For m Xs < 1.1 GeV/c 2 , exclusive B — > K^£ + £~ decays 
are generated with a b — * s£ + £~ decay model according 
to 0,111 . The remaining decays, for mx s > 1.1 GeV/c 2 , 
are generated according to a quark- level calculation |9j 
and the 6-quark Fermi motion model of [Hj- JET- 
SET 0| is then used to hadronize the system consisting 
of a strange quark and a spectator quark. 

The full reconstruction method exploits the 
strong kinematic discrimination provided by 

niES = V E beam ~ Pb and AE = E B ~ E beam , 



where Ebeam is the beam energy and Eb (pb) is the re- 
constructed B-meson energy (three- momentum). These 
quantities are evaluated in the e + e~ center-of-mass 
(CM) frame. 

Clean identification of the B decay products is impor- 
tant for minimizing the backgrounds. Electron candi- 
dates are required to have a laboratory-frame momen- 
tum pi > 0.5 GeV/c and are identified with measure- 
ments from the tracking systems and the electromagnetic 
calorimeter. Bremsstrahlung photons are recovered by 
combining an electron with up to three photons within a 
small angular region around the electron direction |l2| . 
Muon candidates are required to have pg > 1.0 GeV/c 
and are identified as penetrating charged particles in the 
instrumented flux return. Charged kaon identification 
relies on measurements performed with the Chercnkov 
ring-imaging detector and the tracking systems. The 
K° candidates are required to have \m^+^~ — mjfo| < 
11.2 MeV/c 2 , a decay length greater than 2 mm, and 
cos 8 > 0.99, where 5 is the angle between the K® momen- 
tum vector and a line that connects the primary vertex 
with the K® vertex. Charged tracks that do not satisfy 
tight or K ± identification are considered to be pions. 
Neutral pions are required to have a laboratory-frame 
energy greater than 400 MeV, a photon daughter energy 
greater than 50 MeV, and |m 77 — m^ol < 10 MeV/c 2 . 

The B candidates are reconstructed by first selecting 
the e + e~ or pair with the largest value of \p£+\ + 

\pi-\- Then, B — > X s £ + £~ candidates are formed by 
adding any of the following hadronic topologies: K + , 

K+TT°, K+TT-, K+IT-TT , K+ ^ TT+ , K°, K°TT° , ^°7T+, 

KgTT + Tr , and K®ir + ir~. (We name one member of a 
pair of charge-conjugate states to refer to both, unless 
we specify otherwise.) A limit is imposed on the number 
of pions because the expected signal-to-background ratio 
drops significantly with increasing multiplicity. 

With loose selection criteria, the average number of 
B candidates per event is five in the signal simulation. 
A likelihood function for the signal is constructed based 
on the simulated distributions of AE, \og(Psvtx), and 
cos 9b, where Psvtx is the fit probability for the B vertex 
constructed from charged daughter particles and 9b is 
the angle between %>b and the e~ beam axis. We select 
the candidate with the largest signal likelihood before 
further cuts are applied. This approach entails a loss of 
8% in signal efficiency but reduces the overall background 
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by 34% in the final sample. 

Combinatorial backgrounds from nonsignal BB and 
continuum e + e~ — > qq (with 5 = u, d, s, c) events are 
reduced by requiring mx s < 1.8 GeV/c 2 , 5.20 < to_es < 
5.29 GeV/c 2 , and -0.2 < AE < 0.1 GeV. We impose 
a limit on mx s because the signal-to-background ratio 
drops as mx s increases. 

The dominant background producing a peak in m es at 
the B mass originates from B — > J/ipX, ip(2S)X decays 
with J/V> -> £+£" and ^(25) -> £+^". This charmo- 
nium background is suppressed by vetoing B candidates 
with dilepton mass in the ranges 2.70 < m e+e - < 3.25, 
2.80 < to m + m - < 3.20, 3.45 < m e+e - < 3.80, and 
3.55 < m^+p- < 3.80 GeV/c 2 . In the electron chan- 
nel, the veto is applied before and after bremsstrahlung 
recovery to allow for imperfect recovery. The potential 
peaking background from B —> X s ^ decays with conver- 
sion of the photon into an e + e~ pair in the detector ma- 
terial is removed by requiring m e + e - > 0.2 GeV/c 2 . This 
cut also serves to reduce contributions from b —> s e + e~ 
transitions in which the e + e~ pair originates from a pho- 
ton. 

The final suppression of the combinatorial background 
is achieved with a likelihood based on nine variables: (i) 
AE, (ii) AE ROE , (iii) m^ E , where ROE refers to the 
rest of the event (all charged tracks and photon candi- 
dates not included in the B candidate), (iv) the separa- 
tion Az between the two leptons along the beam direction 
measured at their points of closest approach to the beam 
axis, (v) log{PBvtx), (vi) cos0 miss , where 9 mlss is the an- 
gle between the missing momentum vector for the whole 
event and the z axis in the CM frame, (vii) cos#b, (viii) 
cos 9t I , where 8t is the angle between the thrust axes of 
the B candidate and the ROE in the CM frame, and (ix) 
the ratio R2 of the second and zeroth-order Fox- Wolfram 
moments HH- The variables AE, AE ROE , and m§^ E 
are most effective at rejecting BB background, especially 
for events with two semileptonic decays that have large 
missing energy. The event-shape variables | cos 9t | and 
R2 are most effective at suppressing continuum events. 
A likelihood value is computed as the product of nine 
independent probability density functions (PDF) for the 
signal, BB, and continuum background components. Us- 
ing simulated B —* X s £ + £~ decays, we choose cuts on the 
ratio C R = C sianal I (C siBnal + C bB + C cont ) to maximize 
the statistical significance of the signal. This optimiza- 
tion is performed separately for electron and muon chan- 
nels in the regions mx s < 0.6, 0.6 < raj s < 1.1, and 
1.1 < nix, < 1.8 GeV/c 2 , and results in progressively 
harder Cr cuts for increasing mx s ■ 

After applying all selection criteria, we obtain a sam- 
ple of 349 (222) events in the electron (muon) channel. 
According to the simulation, the remaining background 
consists mostly of BB events. 

We perform an extended unbinned maximum- 
likelihood fit to the rriES distribution to extract the sig- 




FIG. 1: Distributions of tties for selected (a) B — ^ X 3 £ + £~ 
(£ = e, n) and (b) B — > X s e ± fT F candidates. The solid lines 
represent the result of the fits and the dashed line represents 
the background component under the signal peak. 



nal yield N S i g as well as the combinatorial background 
shape and yield. The likelihood function consists of four 
components: signal, charmonium peaking background, 
hadronic peaking background (see below), and combina- 
torial background. 

The shapes of the signal and charmonium peaking 
background are described by the same Gaussian PDF, 
with a width of 2.80, 2.61, and 2.74 MeV/c 2 in the 
electron, muon, and electron+muon channels, respec- 
tively. The Gaussian mean and width are determined 
from signal-like B — > J/tpX, ip(2S)X decay candidates 
satisfying all selection criteria but failing the charmonium 
veto (charmonium- veto sample). The charmonium peak- 
ing background is estimated from the simulation to con- 
tribute 0.4 ±0.2 (1.4 ±0.5) events in the electron (muon) 
channel. The size and shape of the hadronic peaking 
background component arising from B — > D^*'nn (n > 0) 
decays with misidentification of two charged pions as lep- 
tons, are derived directly from data by performing the 
analysis without the lepton identification requirements. 
Taking the ir — > £ misidentification rates into account, we 
estimate the hadronic peaking background to be < 0.1 
and 2.4 ± 0.8 events in the electron and muon channels, 
respectively. The error on the misidentification rates de- 
rived from data control samples dominates the uncer- 
tainty. The PDF for the combinatorial background aris- 
ing from continuum and BB events is an ARGUS func- 
tion with an endpoint equal to Eb eam - 

The fit results are presented in Fig.^a) and Table [I] 
The statistical significance is S = \J 2 \n(C ma x / C® nax ) , 
where C ma x is the maximum likelihood for the fit and 
£-max 1S t na t f° r a n t with signal yield fixed to zero. The 
B — > X s £ + £~ signal yield is obtained from a fit to the 
combined electron and muon data. Figure^b) shows the 
ities distribution of B — > X s e ± /i T candidates selected 
using the nominal criteria but requiring leptons of dif- 
ferent flavor. An ARGUS function fits this distribution 
well, consistent with that sample being pure combinato- 
rial background. 
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TABLE I: Signal yield, significance, efficiency, and branch- 
ing fraction in the electron and muon channels, as well as for 
the electron+muon average (first, second, and third rows, re- 
spectively). For the signal yield, the first error is statistical 
and the second error is systematic. For the signal efficiency, 
the first error corresponds to the experimental systematic un- 
certainty arising from detector modeling, hadronization, and 
Monte Carlo statistics, whereas the second error corresponds 
to the uncertainties in the signal model. For the branch- 
ing fraction, the errors correspond to statistical, experimental 
systematic, and signal model systematic uncertainties. 



(%) 



B (1(T 



29.2 ± 8.3 ± 1.3 4.0 2.74 ± 0.27 ± 0.49 6.0 ± 1.7 ± 0.7 ± 1.1 
11.2 ± 6.2 ± 0.9 2.0 1.26 ±0.12 ± 0.25 5.0 ± 2.8 ± 0.6 ± 1.0 
40.1 ± 10.4 ± 1.7 4.3 2.00 ± 0.19 ± 0.37 5.6 ± 1.5 ± 0.6 ± 1.1 
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FIG. 2: Differential branching fraction as a function of (a) 
hadronic mass and (b) dilepton mass, averaged over electron 
and muon channels for data (points) and signal Monte Carlo 
(histogram). The outer (inner) error bars correspond to the 
total (statistical) uncertainties. 



The branching fraction is B = N S i g /(2N BB - e), where 



N bB = (88.9 ±1.0) x 10 6 is the number of BB pairs 
and e is the signal efficiency. The yield N S i g includes a 
contribution from events containing a true B — > X a £ + £~ 
decay but having a selected B candidate with incom- 
plete or wrong decay products (cross- feed events). We 
estimate the cross-feed contribution to the signal yield 
to be 1.5 ± 1.5 (0.6 ± 0.6) events in the electron (muon) 
channel by including in the fit a component whose shape 
and normalization are taken from simulation. The corre- 
sponding uncertainties are estimated with an ensemble of 
simulated data-sized experiments. The signal efficiency 
is adjusted to reflect this contribution. Figure shows 
the B — > X s £ + £~ differential branching fraction as a 
function of mx, and m^+i-, obtained by applying the 
nominal likelihood fit procedure to the data in bins of 
mx 3 and m^+g- . We use the nominal Gaussian shape for 
all bins, as we found no significant shape dependence on 
mx B or m i+e - . 

We evaluate the systematic uncertainties in the sig- 



nal yield by varying the signal Gaussian parameters 
(mean and width) and the signal shape (using asymmet- 
ric signal shapes) within the constraints allowed by the 
charmonium-veto sample. The amount of charmonium 
and hadronic peaking background is varied, as well as 
the shape of the latter. 

Uncertainties in the signal efficiency originate from the 
detector modeling and from the simulation of signal de- 
cays. From control data samples we find uncertainties of 
1.3% per track (0.8% per 7r + ) in the tracking efficiency; 
0.85% per electron, 1.55% per muon, 1.0% per K + , and 
1.4% per 7r + in the charged-particle identification effi- 
ciency; 1.5% per K® and 5.6% per tt° in the reconstruc- 
tion efficiency. We check the efficiency of the cut in Cr 
with the charmonium-veto sample and take the discrep- 
ancy with the simulation (3.7%) as the uncertainty. The 
fraction of signal cross feed included in the signal Gaus- 
sian is varied by ±100%. Altogether, the uncertainty 
from detector modeling is 8.3%. 

The dominant source of uncertainty (18%) arises from 
the decay model. The fractions of B — > K£ + £~ and 
B —> K*£ + £~ decays are both increased (or decreased) 
together by the theoretical uncertainties |l(, resulting in 
the largest contribution (16%). Parameters of the Fermi 
motion model are varied within limits allowed by mea- 
surements of hadronic moments in semileptonic B de- 
cays 1 151 and the photon spectrum in inclusive B — + X s 7 
decays [lfif . The transition point (in mx s ) between the 
B — > K*l + £~ and b — > s £ + £~ decay models is varied by 
±0.1 GeV/c 2 . 

Hadronization uncertainties affecting the region 
mx s > 1.1 GeV/c 2 total 4.9% and are evaluated as fol- 
lows. The ratio between the generator yield for decay 
modes containing a K® and that for modes containing a 
charged kaon is varied within the range 0.50 ± 0.05, to 
allow for isospin violation. Similarly, the ratio between 
modes with one and no ir° is varied within the range 
1.0±0.5, and the ratio between two-body and three-body 
hadronic systems is varied within the range 0.5±0.3. Un- 
certainties in the last two ratios are set by the level of 
dis crep ancy between data and simulation as measured 
in |17| . The uncertainty in the fraction of modes with 
pion or kaon multiplicities different from those used in 
this analysis, or with photons that do not originate from 
7T° decays but rather from 77, 77', cj, etc, is estimated by 
varying these different fractions by ±50%. 

Table [I] summarizes the results of the analysis. For 
the combined B —> X s £ + £~ branching fraction we as- 
sume equal branching fractions in the electron and muon 
channels, and average over both channels. Table IIII shows 
the partial branching fractions in several dilepton and 
hadronic mass ranges. 

We search for CP violation in the B — > X s £ + £~ decay 
by performing separate fits to B and B final states, where 
the final state flavor is determined by the kaon and pion 
charges (modes with X s — K®, K®tt°, or ifg7r + 7T - are not 
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TABLE II: Partial branching fractions in bins of dilepton and 
hadronic mass averaged over electrons and muons. The first 
error is statistical and the second error is systematic. 



m e+e - 
GeV/c 2 ) 



B 

(icr 6 ) 



(GeV/c 2 ) 



B 

(io- 6 ) 



0.2 
1.0 - 2.0 

2.0 - mj/^ 



L0 0.08 ± 0.36+^04 0.4 



mj/4. 



m 



1.6 ±0.6 ±0.5 
1.8 ±0.8 ±0.4 
1.0 ±0.8 ±0.2 



5.0 0.64 ±0.32 



+0.12 



0.6 0.53 ±0.17 ±0.04 

0.6 - 0.8 0.32 ±0.24 ±0.04 

0.8 - 1.0 0.64 ±0.40 ±0.08 
1.0-1.8 3.5 ± 2.1±o;| 



1.0-2.45 1.8 ±0.7 ±0.5 
3.8 - 5.0 0.50 ± 0.25+nnf 



used). We find Nf ig = 14.7±6.5(stat) and N? ig = 22.9 ± 
7.4(stat), corresponding to an asymmetry Aqp = (N B — 
N B )/(N E + N B ) = -0.22 ± 0.26(stat). The systematic 
uncertainty is taken to be the statistical uncertainty in 
the asymmetry measured in the charmonium-veto sample 
Agp = -0.005 ± 0.016(stat). 

In summary we observe a signal of 40 ± 
10(stat) ± 2(syst) B -> X s £+£- events with 
a statistical significance of 4.3 a. The corre- 
sponding branching fraction B(B — > X s £ + £~) = 
(5.6 ± 1.5(stat) ± 0.6(exp syst) ± l.l(model syst)) x 
10 -6 for m£+f- > 0.2 GeV/c 2 agrees well with predic- 
tions and a previous measurement In restricted 
dilepton mass ranges below and above the charmo- 
nium region (see bottom of Table [H}, the partial 
branching fractions also agree with predictions (lif . 
At low hadronic mass, the partial branching fractions 
are consistent with existing B — > K^£ + £^ measure- 
ments 0, H| . We determine the direct CP asymmetry to 
be A C p = -0.22 ± 0.26(stat) ± 0.02(syst), in agreement 
with a p redicted asymmetry of (0.19!q '{I) x 10~ 2 in the 
SM 191. Overall, we find no evidence for contributions 
from physics beyond the SM. 
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